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Abstract 

We consider a model of restricted dimers coupled to two-dimensional causal dynamical tri- 
angulations (CDT), where the dimer configurations are restricted in the sense that they do 
not include dimers in regions of high curvature. It is shown how the model can be solved 
analytically using bijections with decorated trees. At a negative critical value for the dimer 
fugacity the model undergoes a phase transition at which the critical exponent associated to 
the geometry changes. This represents the first account of an analytical study of a matter 
model with two-dimensional interactions coupled to CDT. 
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1 Introduction 



The program of Dynamical Triangulations (DT) was an attempt to make the notion of 
quantum gravity as the path integral over all Euclidean geometries well defined (see [1] for 
reviews). The idea of DT is to define the path integral by a sum over discrete approxima- 
tions to each geometry. The approximations are formed by gluing together basic polygonal 
building blocks; such blocks are often taken to be triangles, or higher dimensional analogues, 
and therefore the discrete geometries are referred to as triangulations. Once the sum over 
the discrete structures is computed, a scaling limit must be taken in which the mean number 
of triangles in a discretisation diverges while simultaneously allowing the size of the triangles 
to go to zero. Intuitively such a limit should return one to the case of a continuous geometry. 

In two dimensions this program was quite successful for a number of reasons. Most 
importantly, for two dimensions the dual graph to the triangulations can be generated as 
the Feynman diagrams of matrix models and hence the sum over such graphs is reduced to 
the problem of computing a certain matrix integral, for which many techniques are known 
(see for instance P). Furthermore, the matrix model approach allowed for the counting of 
decorated triangulations which could be interpreted as statistical systems such as the Ising 
model coupled to the quantum geometry [2]. Another reason for the original interest in two- 
dimensional DT is that there exists a purely continuum approach to the two-dimensional 
gravitational path integral, in which its computation is reduced to the study of Liouville 
theory [31 S] . There was thus a way of directly comparing the scaling limit of DT with the 
continuum calculation and in all greement was found. In particular the decorated 

graphs were found to scale to {p, q) minimal model conformal field theories (CFTs) coupled 
to Liouville theory |3l HI |5]. Because all such CFTs have a central charge satisfying, c < 1, 
the inability to analytically analyse DT coupled to c > 1 matter became known as the c = 1 
barrier. Numerical simulations suggested the reason for this barrier was that the geometry 
of the space degenerated into a branched polymer phase beyond this limit [6]. 

In higher dimensions, the program of DT was less successful (we refer to [Tj for an 
overview and references): Depending on the value for the gravitational constant two distinct 
phases were found; the crumpled phase, with a Hausdorff dimension du ~ oo, and the 
branched polymer phase, with Hausdorff dimension dn = 2. Furthermore, the transition 
between these phases was found to be of first order and therefore there was no possibility 
of finding a critical point at which an extended four- dimensional geometry would emerge. 
To rectify this, a proposal was made in [7j to instead consider the sum over a restricted 
class of triangulations, known as causal triangulations, in which a preferred time foliation 
was present and spatial topology change was forbidden with respect to this slicing. This 
approach is known as Causal Dynamical Triangulations - CDT for short. 

An unfortunate aspect of the casual restriction of CDT was that the relation to matrix 
models, as well as to a continuum path integral formulation was lost - although there has 
been progress in establishing a matrix model formulation [8l [9l |T0] as well as relations to a 
continuum path integral in the proper time gauge [H] (see also [12] for a review). Due to the 
loss of the matrix model formulation, the original study used a transfer matrix method to 
solve the pure case in two dimensions and in higher dimensions all work has been numerical 
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in nature (see [13^ for a recent review). This has the consequence that unhke DT, it has 
been very difficuh to analytically investigate matter coupled to CDT. One exception to 
this was the work of [H], in which a peculiar form of dimer matter was studied in which 
dimers were allowed to be of arbitrary length in the timelike direction. Although they were 
able to obtain a multicritical point, this model was in some sense one- dimensional, due to 
the absence of spacelike interactions between dimers and it is unclear what sort of theory 
the continuum limit represents. Despite this absence of analytical work, there have been 
a number of very interesting numerical studies of the behaviour of various forms of matter 
coupled to two-dimensional CDT [15]. 

These studies indicated two intriguing properties; firstly, unlike DT, CDT does not appear 
to modify the critical exponents of the CFT and secondly, the c = 1 barrier appears to be 
of a milder form with indications the geometry instead changes to a "blob" type geometry 
of higher dimension [TBj|. These properties are currently conjectures based on numerical 
evidence and therefore it is a pressing open problem in the study of CDT to find a method 
by which matter coupled to CDT can be analysed analytically. 

In this paper we introduce a method of analytically computing the partition function for 
a certain kind of hard dimer model coupled to CDT. Our method is in the spirit of the tree 
bijection techniques used by [17] in which it was shown that conformal matter coupled to 
DT could be analysed using such bijections without ever referring to matrix integrals. 



The particular matter model we will consider is known as hard dimers. For a regular tri- 
angular lattice a hard dimer model consists of objects known as dimers which occupy two 
adjacent sites in the lattice. Each dimer present in the lattice contributes a factor of ^ to the 
weight associated to the entire configuration. The interaction between the dimers is via a 
"hard" interaction in which only a single dimer can reside inside a given triangle. This prop- 
erty makes dimers somewhat easier to analyse than other models such as the Ising model, 
as the interaction only manifests itself through the number of valid configurations. 



To couple dimers to DT we make use of matrix model techniques. Since matrix models 
generate the dual graph of the triangulations it is easiest to proceed by reformulating the 
dimer model on the dual graph. Such a reformulation is simple; the model corresponds to 
colouring the links of the graph with two colours, black and red, such that red links can 
never end on the same vertex. The red links then correspond to the dual of the dimers. 

A model which generates Feynman diagrams appropriate for counting such a system is 
given by the Hermitian two-matrix model for matrices of size N x N, 



2 Hard Dimers in DT and CDT 



2.1 Dimers in DT 




(1) 
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where A generates red links and M black. In this model g is the fugacity associated to a single 
triangle and ^ is the weight associated to the inclusion of a dimer. In the continuum limit g 
will scale to a coupling describing a perturbation by the lowest dimension operator present 
in the continuum theory. Usually such an operator is the identity operator and therefore g 
may be interpreted as the cosmological constant A. However, for non-unitary theories there 
exist negative dimension operators and g can no longer be given this interpretation. 

Proceeding from ([1]) it is clear that we may integrate out A to obtain, 

[rfMje-^^'^^^^^) (2) 



where V{M) is a polynomial of order four. This model was among the first matrix models 
to be analysed in the context of two-dimensional quantum gravity [18j where it was found 
to contain a critical point {gc, ^c) about which a scaling limit could be defined. Initially, 
the scaled theory was thought to correspond to the Ising model coupled to gravity, however, 
it was later realised [H] that it actually corresponded to the critical point of the hard 
dimer model, described by the (2, 5) minimal CFT. The reason for this confusion was that 
the original work [18] used ([2]) rather than ([1]), thereby obscuring the interpretation of the 
degrees of freedom and so the identification of the scaling limit with a particular model was 
done on the basis of matching the observables. In particular, the observables of this model 
are 7, known as the string susceptibility, and defined as the first non-analytic power in, 

Z ~ A^-^, (3) 

and the critical exponent a, defined similarly through ^19] , 

d^oggdO 



- ^cT. (4) 



Note that for models in which the identity operator is the lowest dimension operator, ([3]) 
can be written as, 

Zr^ig-g,y-\ (5) 

For models in which the identity is not the lowest dimension operator we call the exponent 
associated with the modified string susceptibility. The string susceptibility characterises 
the modification of the geometry by the presence of matter. It was this observable that was 
used to identify the scaling limit of ([2]) about its multicritcal point with the Ising model, due 
to false identification of the string susceptibility with its modified version, which happened 
to have the same value in this case. By solving the model ([2]) it was found that, 

7 = 7modified = --, ^ ^ 2' 

This should be contrasted with the case of dimers on a fiat regular lattice, for which cr = — |. 
The difference is described by the famous KPZ formula [1]. 



4 



Figure 1: A configuration of a CDT coupled to dimers as considered in [20]: Tliere are no 
interactions between adjacent spatial slices. 

2.2 Dimers in CDT 

There has been one previous attempt to study CDT coupled to the same sort of hard dimer 
matter as found in DT [20]. However, the result of this study was that although the discrete 
model could be solved, the scaling limit was still in the universality class of pure CDT. Before 
moving on to our own construction it is instructive to understand the failure of this previous 
study. In [20j a generalisation of CDT was analysed in which a discrete version of a higher 
curvature term was added to the measure. This corresponded to attaching an extra weight 
to each instant in a triangulation when two adjacent triangles in the same space-like slice 
are of the same orientation. It was shown in [20] that this weighting maps to a dimer-like 
system in which the dimers are confined to lie only across time-like edges connecting an up- 
and a down-pointing triangle, resulting in configurations similar to that shown in Figure [H 

For statistical models such as the Ising or dimer model, the critical point, at which the 
model is described by a conformal field theory, corresponds to the situation in which the 
correlation length diverges. It is therefore perhaps not surprising that the dimer-like model 
considered in [20] was still in the universality class of pure gravity; there is no correlation of 
dimer configurations even between adjacent layers. This effectively transform the problem 
into a succession of independent dimer models on one-dimensional circles, for which it is 
known that no phase transition exists. It is clear then that in order to move out of the 
universality class of pure CDT we must at least allow for some interaction between adjacent 
layers. This is the subject of the next section. 

3 Decorated Trees and Dimers 

It has been known for some time beginning with the work of Cori-Vauquelin [21] and Scha- 
effer [22] that there exist bijections between tree graphs and various kinds of triangulations 
(not necessarily composed of triangles), known in the mathematical literature as maps. A 
particularly interesting use of these bijections have been the solutions of the Ising model and 



5 



Figure 2: A causal triangulation drawn 
in the plane. The dashed links emanat- 
ing from the outermost circle all end 
on a single vertex in the infinite face, 
which we have omitted from the dia- 
gram. 



Figure 3: Casual triangulations have a 
bijection with trees. In this case the 
tree associated to the triangulation is 
shown by thick lines. We have also 
omitted the root vertex attached to the 
central vertex. 



other models [17] coupled to two-dimensional gravity without the use of matrix models. The 
purpose of this paper is to apply these techniques to the problem of hard dimers coupled to 
CDT. 



3.1 A Tree Bijection for Pure CDT 

The bijections formulated in [22| [23] related maps to various kinds of labelled trees. The 
maps themselves were precisely of the kind that contribute to the partition function of DT. 
To apply these techniques we first review a tree bijection from CDT to rooted planai0 trees. 
This bijection was introduced in [20] and has been used for example to bound the spectral 
dimension of the spacetimes resulting from CDT [24j. The bijection is most easily described 
by first fiattening the CDT graph for a cylinder on to the plane, as is shown in Figure [2j 
Next, we denote the vertices in the kth spacelike slice as Sk, we define 5*0 as containing a 
single new vertex connected to all vertices in 5*1. The following steps then map the CDT to 
a rooted planar tree T; 

1. All vertices in the triangulation are in T in addition to a new root vertex, called the 
root, that is only linked to 5*0. 

2. All links from 5*0 to 5*1 are in T. 

^Here "planar" refers to the fact that the embedding of the tree in the plane is significant, in that it 
determines an order on the vertices in that one proceeds another in an (anti) clock-wise direction in the 
plane. 
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Figure 4: The dimers configurations 
that are permitted in the restricted 
dimer modeL 



Figure 5: The dimer configurations 
that are excluded in the restricted 
dimer model. 



3. All links from a vertex in Sk to vertices in Sk+i, apart from the anticlockwise-most 
link, are in T. 



These steps produce the tree shown in Figure [31 Denoting the vertices at a height n above 
the root by Vn, this procedure can also be shown to be invertible by the following steps, 



1. Delete the root node and the link joining it to Vi. We identify Vi as Sq. 

2. For all n < h, where h is the total height of the tree, join each vertex in Vn+i to the 
next vertex in the clockwise direction, as determined by the embedding in the plane. 

3. For each vertex v in Vn, 1 < n < h, that is not of order one in the tree, add a link 
between v and a vertex in l^+i such that the new link is the most anticlockwise link 
emerging from v and furthermore maintains the planarity of the graph. 

4. For vertices of order one, add a new link between it and the unique vertex in the next 
layer such that the planarity of the graph is maintained. 



3.2 A Tree Bijection for a Restricted Dimer Model 

We now present an extension of the above bijection to allow certain dimer configurations to 
be encoded in the tree. Given a triangulation, we first label all vertices with the number zero. 
Then consider Figure HI it displays a number of possible ways a dimer could be embedded 
into a triangulation. Note that the list of how a dimer may be embedded in the triangulation 
is not exhaustive and neglects the two possibilities shown in Figure |5l We refer to the model 
which neglects these possibilities as the restricted dimer model. Note that our model includes 
the dimer configurations of [20] as a subset in addition to allowing dimers to lie between 
layers, as is required to move out of the universality class of CDT. 

When one of the restricted dimer configurations appears in a triangulation we will label 
the vertex shown in Figure H] with a number corresponding to how the dimer is configured. 
By noting that the tree bijection of the previous section does not delete any vertices, we can 
extend it to the present case simply by defining it to preserve the labelling of the vertices, 
thereby obtaining a labelled tree. Figure |6l 
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Figure 6: An example of the labelled 
tree encoding the dimer configuration 
in the CDT. 



Figure 7: Here, in yellow, we see two 
invalid labels in the tree. At the top of 
the diagram we see why a type 2 label 
may not have a type 1 label as its anti- 
clockwise most child. At the bottom 
of the diagram we see how the label 3 
causes there to be non-local constraints 
on the labelling. 



We now address the question of what class of labelled trees is produced by the above 
map. Firstly, not every labelled tree is in the range of the map. This fact can be seen by 
noting that certain labelled trees can not arise as they would correspond to a configuration 
of dimers in which some occupied the same site, as shown in Figure [7] by the yellow dimers. 
We therefore have a tree for which certain restrictions are present on the labelling. The task 
is now to characterise these restrictions and furthermore find a way to count the number 
of trees appearing in this class. Note that m doing so we have to account for the non-local 
interactions between different branches of the tree. 

We begin by considering the labelling on a tree. We see that a vertex labelled with a 0, 
1 or 3, places no restrictions on the labelling of any of its children. However, the label 3 will 
affect the next branch in an anti-clockwise direction as can be seen with the lower yellow 
dimer in Figure [71 This effect is non-local in the tree graph and we postpone a discussion on 
how to handle it until later. For a vertex carrying the label 2, we see from the upper yellow 
dimer in Figure [7] that it must have at least one child and furthermore this child can only 
be labelled with a or 2. Furthermore such a dimer does not place any other constraints 
on the labelling of the graph and hence its effect is local in the tree graph. Since the rules 
associated with labels 0, 1 and 2 only put restrictions on the child nodes of a vertex in the 
tree, trees containing only type 0, 1 and 2 dimers can be counted recursively via generating 
functions. 

The type-3 dimers may be included by mapping them onto type-2 dimers via another 




Figure 8: The move used to transform a 
type 2 dimer into a type 3 by changing 
the triangulation. This move can be 
performed independently of the rest of 
the triangulation. 




Figure 9: A graphical representation of 
the recursion relation needed to com- 
pute the generating function for the la- 
belled trees. 



bijection. The type-2 and type-3 dimers have the special property that by flipping the central 
diagonal we may convert one in to another, as in Figure [H and crucially, this does not affect 
any other dimers in the configuration. To clarify, by performing the move in Figure [8] we map 
one triangulation with a particular configuration of dimers to a distinct triangulation with a 
distinct configuration of dimers containing one less type-2 dimer and one more type-3 dimer. 
Hence, by counting each type-2 dimer twice we are able to account for all configurations 
including type-3 dimers. This is a crucial non-trivial observation necessary to perform the 
calculation in this paper |§ 

We are now in a position to count the number of labelled trees satisfying the above rules. 
Let us introduce the generating functions Gi defined by, 

oo oo 
m=0 n=0 

where J\fi{n, m) is the number of trees with a root carrying the label i, containing n vertices, 
m of which carry a non type-0 label, thereby indicating the presence of m dimers. We can 
obtain an expression for Go of the form, 

Go{^, g) = . . = 1 + /(Go + + 2^0^) + ..., (8) 

which is shown pictorially in Figure [91 Given that the type-1 label places no restriction on 
the children we also have Gi = Gq. Lastly, by recalling the restrictions due to a type-2 
dimer, we can write, 

G2i^,9) = g'iGo + 2^G2)Go. (9) 

Solving for Go results in the equation, 

2^?"^'^-^ - g\l + 30GI + (1 + 2g'OGo -1=0. (10) 

The critical line gcriO of this model, at which the mean size of the triangulations diverges, 
can be found from the radius of convergence of Gq. This is easily found by computing the 
discriminant of ( fTOl) . obtaining the following equation for gcriO^ 

GAegl - 4e'(l + 24e + I2e)gi + 4(1 + 8^ + 12^' + 3e)9l - 1 - 6^ - = 0. (11) 

^For the sceptical reader, we have verified that our generating function derived via this argument repro- 
duces the correct number of restricted dimer configurations for triangulations of up to size 18 as found by 
computer. 
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Figure 10: A plot of the value of gcr- Its real value is shown in blue and gold, corresponding 
to the two branches identified at ^ = 0. The imaginary value of the respective branches is 
shown in red and green. One can see that the critical line ends at the value of ^cr = ~(27)^2, 
marked by the dot. 

From the discriminant of the above equation there exist two possible critical values of ^; 
^cr = ±(27)~^/^. Furthermore we identify the correct branch of gcr by requiring that it 
reproduces for ,^ = the known critical value of 5'cr(0) = 1/2 for pure CDT. 

By considering how the roots of ( fTTl) change as ^ varies, one can see that there exists 
one branch that approaches (7cr = 1/2 as —t- 0"^ and another which approaches Qcr = 1/2 
as — )► 0^, with each solution diverging and becoming complex upon passing through = 0. 
If we identify these two branches at ,^ = then we obtain the plot of Qcr show in Figure [TOl 

Expanding the solution in a taylor series around its endpoint at C,cr we obtain using (jlj), 

. = i (12) 

and furthermore substituting the critical value ^ into the solution for Gq we obtain, 

7modificd ~ 2' (-'-3) 

from which it is clear that we are now outside the universality class of pure CDT! These 
results are summarised in Table [1] alongside their counterparts in DT. These results are quite 
striking when compared to the flat space and DT exponents. Firstly, the value of 7modified is 
positive. For all minimal models coupled to DT 7 is negative and historically a positive 7 
was associated with going beyond the c = 1 barrier [25] . 

Secondly, a has exactly the same value as obtained in DT. This is particularly surprising 
given the numerical evidence that coupling to CDT preserves the flat space exponents. One 
should note however that all numerical work on matter coupled to CDT has used matter for 
which c > 0. Of course, one could argue that we have not solved the full dimer model and 
that we might obtain agreement with the numerical work once we include all configurations. 
However, we will now argue that even in the unrestricted case we will have a differing from 
its fiat space exponent. 
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a 


7 


7modified 


Dimers with DT 


1 

2 


3 
2 


1 

^ 


Restricted Dimers with CDT 


1 
2 


? 




Pure CDT 




1 
2 





Table 1: A comparison of the critical exponents of the restricted dimer model coupled to CDT 
and the normal dimer model coupled to DT. In both cases we are outside the universality 
class of the pure theory. 

Note that the dimers shown in Figure [5] correspond to dimers placed in regions of high 
curvature. One consequence of this is that highly curved spacetimes will be entropically 
suppressed. Furthermore, the negative fugacity of the dimers at their critical point will also 
enhance the weight associated to spacetimes that can support dimers. Hence there is a 
double enhancement of the flatter spacetimes. Hence, if our results for the critical exponents 
were to change for the unrestricted dimer model, we would expect them to move further 
from their fiat space counter parts. Therefore, we expect that the full dimer model coupled 
to CDT belongs to the same universality class as the restricted dimer model coupled to CDT 
as considered here. The results are summarised in Table [H along side their counterpart in 
DT. 



4 Discussion 

In this paper we have solved analytically a hard dimer model with two dimensional in- 
teractions coupled to CDT. We have been able to compute the disc function, the string 
susceptibility and the critical exponent for the hard dimer model, a. Both exponents were 
seen to differ from the pure CDT case, hence indicating we have been able to analytically 
study a known matter model coupled to CDT for the first time. Furthermore, the value of a 
we found matched that obtained from DT. This is particularly surprising given the numerical 
evidence that matter coupled to CDT does not have critical exponents that differ from their 
fiat space values. 

The particular hard dimer model we considered is an extension of the one considered 
in [20] with the crucial difference that we are able to include dimers laying both within a 
single spacelike slice and those which lie across two. The ability to include both types of 
dimers is the fundamental reason that this model lies in a different universality class from 
pure gravity; dimers can interfere with one another across layers, allowing for long range 
correlations to arise. 

We build on the work of |1^ in which bijections were constructed between labelled 
tree graphs and DT with matter, in addition to [20] which introduced a bijection between 
unlabelled trees and pure CDT. In our case we demonstrated that in the case of CDT we 
could keep track of the dimers by introducing a labelling on the tree. This in itself is not 
surprising, however since the connectivity of the tree graph is lower than the triangulation, 
interactions between neighbouring vertices in the triangulation generically become non-local 
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interactions in the tree. A novel aspect of our work is that we were able to account for the 
non-local interactions appearing in the restricted dimer model by introducing a mapping 
between trees. This map relates trees with only local interactions to trees with non-local 
ones, allowing us to count the latter. 

The dimer configurations that we have not included in this model correspond to dimers 
embedded in regions of high curvature, thereby suppressing spacetimes which are highly 
curved. We therefore physically expect our model to scale to a continuum model in which 
there exists a coupling between matter and curvature. Furthermore, we have argued that 
this implies that the value of a will not return to its flat space value upon the inclusion of 
all dimer configurations. 

Note 

While completing this article we were informed by Jan Ambj0rn about an unpublished work 
of him, L. Glaser, A. Gorlich and Y. Sato regarding the multi-critical analogue of a new 
continuum limit of matrix models as introduced in [9] in which they obtain the same critical 
exponents as those in this article. This is very interesting, as the approach of their paper is in 
some sense highly complementary to ours. Indeed their work provides very strong evidence 
that the restricted dimer model belongs to the same universality class as the full dimer 
model, as we argued in this paper. An obvious advantage of the matrix model approach is 
that it is much easier to generalise to other multicritical points, however even in the case 
of DT it was difficult to be certain what continuum theory a particular multicritical point 
scaled to. The identification of a particular scaling limit with a CFT coupled to gravity 
was performed by matching the exponents with those from the KPZ formula. However in 
CDT no such formula exists and therefore our work complements the work of Ambj0rn et al. 
nicely in that it gives an explicit demonstration that the multicritical point found is indeed 
the continuum theory associated with the hard dimer phase transition. 
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